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Bi-Hamiltonian nature of the equation 

t^tx ^xy t^y ^yy ^x 
V. Ovsienko 

Abstract 

We study non-linear integrable partial differential equations naturally arising as bi-Hamiltonian Euler 
equations related to the looped cotangent Virasoro algebra. This infinite-dimensional Lie algebra (con- 
[T | . structed in [16 J is a generalization of the classical Virasoro algebra to the case of two space variables. 

| Two main examples of integrable equations we obtain are quite well known. We show that the rela- 

tion between these two equations is similar to that between the Korteweg-de Vries and Camassa-Holm 
equations. 
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1 Introduction 

The differential equation 

00 ■ 

where u = u(t,x,y) and where u x ,u y , etc. are the partial derivatives, is a nice example of a 
| non-linear integrable model. This equation is quite well known and appears in the Martinez 

' Alonzo-Shabat "universal hierarchy" (see [13] . formula (8)). 

The main purpose of this note is to show that this equation (coupled together with another 
differential equation, see formula (13.61) below) naturally appears as a bi-Hamiltonian vector field 
(in particular, the variable t plays the role of time while x, y are space variables). More precisely, 
^ we will show that this equation is an Euler equation on the space dual to the "looped cotangent 

Virasoro algebra" introduced in [16] . This, in particular, implies its integrability in the (weak 
algebraic) sense of existence of a hierarchy of first integrals in involution. 

The bi-Hamiltonian approach to the same Lie algebra has already been considered in [16] 
and led to another non-linear differential equation: 

ft = fx d^fy - fyU + Cd' 1 fyy, 

which can be rewritten without non-local terms: 

Utx — l^xx l^y t^xy ^x CLlyyt (1*^) 

after the substitution f = u x . Here c S K is an arbitrary constant (the "central charge" ) . Note 
that this equation is also a quite well known integrable system (see [5j [6] and also [1]) that 
appears both in differential geometry and hydrodynamic. 



Equations (|1,2|) and look alike but they are not equivalent to each other. We will show 
that the relation between these equations is similar to that between the classical Korteweg-de 
Vries equation (KdV) and the Camassa-Holm equation (CH). Recall that both KdV and CH 
are bi-Hamiltonian systems on the dual of the Virasoro algebra, see |12|. [31 \TT\ [T4l HOj , An 
interesting "tri-Hamiltonian" viewpoint was suggested in [15], in order to establish a certain 
duality between KdV and CH. Equations (jl.2p and (jl.ip are dual in the same sense. 

This paper fits into the general framework due to V.I. Arnold, see [l]. Non-linear partial 
equations are viewed as Euler equations on the dual of a Lie algebra (for instance, the Lie 
algebra of vector fields). This approach explains the geometric meaning of the equations: every 
Euler equation describes geodesies of some left-invariant metric on the corresponding group (of 
diffeomorphisms) . 

The bi-Hamiltonian Euler equations are of special interest. Most of the known bi-Hamiltonian 
non-linear partial differential equations (KdV,CH, etc.) are of dimension 1 + 1 (i.e., contain only 
one space variable). Equations (jl.ip and (|1.2I) provide with examples of such equations in the 
(2 + l)-dimensional case. 

2 The bi-Hamiltonian formalism on the dual of a Lie algebra 

In this section we recall the general construction of pairs of compatible Poisson structures on 
the space dual to a Lie algebra. We also give the standard construction of bi-Hamiltonian vector 
fields on this space, due to F. Magri [12]. 

Let a be a (finite-dimensional) Lie algebra, the canonical Lie-Poisson(-Berezin-Kirillov- 
Kostant) bracket on a* is given by 

{F,G}{m) = ([d m F,d m G],m), (2.1) 

where m € a* and where d m F and d m G are the differentials of F and G at m understood as 
elements of a, namely dF m 6 (a*)* = a. This Poisson structure is linear, i.e., the space of linear 
functions equipped with the bracket (|2.1j) is a Lie subalgebra of C°°(a*) (isomorphic to a). 

Given a skew-symmetric bilinear form u : a A a — ► R, one defines another Poisson structure 
on a: 

{F,G} u) (m)=uj(d m F,d m G). (2.2) 

This structure is with constant coefficients, i.e., the bracket of two linear functions is a constant 
function on a*. 

Two Poisson structures are called compatible (or a Poisson pair) if their linear combination 
is again a Poisson structure. The following simple fact is well known (see, e.g., [2], Section 5.2). 

Proposition 2.1. The Poisson structures 112.1]) and 112.2]) are compatible if an only if uj is a 
2-cocycle on a. 

The simplest example of a constant Poisson structure (|2.2p corresponds to the case where 
the 2-cocycle to is trivial (i.e., a coboundary). Every such structure is of the following form. Fix 
a point mo £ a* and set 

u)(x, y) = (m , [x,y]) . (2.3) 

It worth noticing that one can understand this particular case of constant Poisson structure on 
a* as the most general one. Indeed, it suffices to replace a by its central extension. 

Every function H on a* defines two vector fields that we denote Xh and Xfj on a* : the first 
one is Hamiltonian with respect to the linear structure (|2.ip and is given by 

X H (m) =ad* dmH m, (2.4) 
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while the vector field is Hamiltonian with respect to the constant bracket (|2.2|) , In the 
particular case (|2.3p . one has explicitly 

X%(m)= a d* dmH m . (2.5) 

Given two compatible Poisson structures, a vector field which is Hamiltonian with respect to 
the both structures is called bi-Hamiltonian. The usual way to construct bi-Hamiltonian vector 
fields on a* is as follows. Consider the following 1-parameter family of Poisson structures 

{,}a = {,}.-A{,}, 

(parameterized by A € R). Assume that H is a Casimir function of this bracket, i.e., one has 

{H, F}\ = 0, for all F € C°°(o*). 
Assume also that H is written in a form of a series 

H = H + XH 1 + X 2 H 2 + --- (2.6) 

One immediately obtains the following facts: 

1. the function Hq is a Casimir function of { , } w ; 

2. the Hamiltonian vector field corresponding to are bi-Hamiltonian, namely 

x H k = X% k+1 , 

for all k; 

3. all the functions Hf, are in involution with respect to the both Poisson structures, indeed, 
for k < i one has 

{Hk, Hi} = {iTfc+i, He} w = {-fffc+i, He-i} = • • • = 0, 

and therefore are first integrals of every vector field Xjj k ■ 

Let us summarize the method. To construct an integrable hierarchy, one chooses a function 
Hq which is a Casimir function of the constant Poisson structure { , } w ; one then considers 
its Hamiltonian vector field, Xjj , with respect to the Lie-Poisson structure. This vector field 
is again Hamiltonian with respect to the constant Poisson structure, with some Hamiltonian 
function H±, so that one has: Xjj = Xjj . One then iterates the procedure to find H 2 , H%, etc. 

3 The looped cotangent Virasoro algebra and its dual 

In this section we recall the definition [16] of the looped cotangent Virasoro algebra. We also 
describe its coadjoint representation. 

Let us start with the definition of the classical Virasoro algebra. Consider the Lie algebra, 
Vect(5 1 ), of vector fields on the circle: f(x) J| where / € C°°(S 1 ) and x is a coordinate on S 1 , 
we assume x ~ x + 2tt. To simplify the formulas, we will identify Vect(S' 1 ) with C°°(S l ); the 
Lie bracket in Vect(5 1 ) is then given by 

[/, g]= f 9x~ fx 9- 
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The Virasoro algebra, Vir, is a (unique up to isomorphism) one-dimensional central extension 
of Vect(S 1 ). It is defined on the space Vect(5 1 ) © R, the commutator being given by 



[(/, «), (g, P)] = fg x - f x g, / f gxxxdx 



(3.1) 



Note that the constants a and (3 do not enter the right hand side of the above formula since 
they belong to the center of Vir. 

The Virasoro algebra was found by Gelfand and Fuchs [8], the constant term in the right 
hand side of (|3,ip is called the Gelfand-Fuchs cocycle. This Lie algebra plays an important 
role in mathematical physics, essentially because of the applications of its representations to 
conformal field theory, but also because of its applications to integrable systems. 

The dual space, Vect(5 1 )*, is the space of distributions. One often considers only a subspace, 
Vect(5 1 )* eg , called the "regular dual" (cf. [9|). As a vector space, this regular dual is, again, 
isomorphic to C°°(5' 1 ), the pairing (., .) : Vect(5 1 ) <g> C°°(S l ) — > R being give by 



d 



f{x) a(x) dx. 



The regular dual to the Virasoro algebra is Vir* eg = C°°(S 1 ) ©1 
its regular dual is: 



the coadjoint action of Vir on 



ad 



(/,«) 



(a, c) = (/ a x + 2 f x a + c f xxx , 0) . 



This formula easily follows from (I3.ip and the definition of ad*, see [9]. Note that the constant 
c is preserved by the action, it is therefore a parameter called the central charge. 

Remark 3.1. The Virasoro algebra is, indeed, exceptional. The reason is that the Lie algebras 
of vector fields on a manifold of dimension > 2, has no central extensions, cf. [7]. The problem 
of generalization of the Virasoro algebra is an interesting subject studied by many authors. 



The looped cotangent Virasoro algebra [16] is a generalization of Vir in the case of two 
variables. We consider the 2-torus T 2 and define a Lie algebra structure on the space 



= C oo (T 2 )ffiC oo (T 2 )©: 



given by the commutator 

// \ (9 

b 

\ (a, a') J \ (f3,(3') ) 



( fgx- fxg 

fb x + 2f x b 



V 



fgx 



ga x -2g x a 

x dxdy, \ 

JS 1 xS l 



(fb y -ga y ) dxdy 



(3.2) 



where (x, y) are the usual coordinates on T 2 and where /, g, a, b are smooth functions in x, y; 
the constants a, a' , /3, G R are elements of the center. Note that, unlike the Virasoro algebra, 
the center of g is two-dimensional. 



Remark 3.2. One notices that the quotient-algebra 

with coefficients in the semidirect sum Vect(5 1 ) x Vect(5 i )* eg . The dependence in y-variable 



2 (by the center) is the loop algebra 



in this quotient-algebra is somehow trivial. The second 2-cocycle in (|3.2|) . however, makes this 
dependence in y non-trivial. Note also that this cocycle is rather similar to the Kac-Moody 
cocycle. 
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We will need the coadjoint representation of g and the notion of regular dual space. Consider 
the pairing (., .) : g ® g — > R 





\ 


(9 \ 


a 




b 


\ (ai,«2) 


/ 


\ (ai,a 2 ) / 



/ (f b + g a) dxdy + ai(3i + 0.202, 
Js 1 xS 1 



that identifies g with a part of its dual space: g g*, we call this subspace the regular dual 
space of g and denote it by g* eg . The coadjoint action of g on g* og can be easily calculated: 



ad 



( f ) 

a 

V («1>«2) / 



(9 \ 




b 


H 


\ (ci,c 2 ) J 





f b x + 2 f x b - a x g - 2 a g x + c\ f xxx + c 2 a y 



(3.3) 



Note that the center R 2 C g acts trivially. 

The Lie algebra g is infinite-dimensional. In order to define the brackets (|2.ip and (|2.2p in 
this case, we consider only the space of so-called pseudodifferential polynomials on g* cg : 

H(f, a) = h (/, a, f x , a x , f y , a y , d" 1 /, d~ l a, d' 1 /, d~ l a, f xy , a xy , . . .) dxdy, 

JS 1 xS 1 

where h is a polynomial and /, a, f x , a x , f y , a y , . . . are understood as independent variables. 

The differential d m H is replaced by the standard variational derivative: 

d (f>a) H:= (5 a H, S f S) 

understood as element of g/M 2 . The Lie-Poisson structure (|2.ip then makes sense on g* eg and 
the Hamiltonian vector fields are again given by (12. 41) . 

Example 3.3. Recall that the Euler-Lagrange equation provides an explicit formula for varia- 
tional derivatives. For instance, one has 

5 a H = h a -8 x (h ax ) - 3y (h ay ) - (hg-l a ) - dy 1 (hg-l a 



+ (d X ) 2 (h axx ) + d x d y (h axy ) + (d y ) 2 (h a 



where, as usual, h u means the partial derivative similarly h ax = etc.. 

One of course should be careful with the definition of the non-local operators 9" 1 and dy 1 . 
We use the expression 



(3*V)(s,y) 



f(Z,y)<% 



2k 



f(x,y)dx, 



and similarly for dy 1 . 

We refer to [3] for further details on Hamiltonian formalism on infinite-dimensional (func- 
tional) Lie algebras. 
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3.1 Calculating the bi-Hamiltonian equations 

Let us fix the following point of Q^ cg : 

mo = {fix), a(x), ci, c 2 ) = (1, 1, 0, c) , (3.4) 

with arbitrary c G R, and consider the constant Poisson structure (|2.2p corresponding to the 
coboundary fl2.3|> . The Hamiltonian vector field with the Hamiltonian H is then given by 

ft = -(5 a H) x + c(5 a H) y 

a t = 2(S a H) x -(S f H) x + c(S f H) y . 

The limit case c — > oo corresponds to the following structure 

/t = 

(3.5) 

= 

We are ready to formulate our main result. 
Theorem 3.4. The following system on g* cg 

Utx — l^xy ^yy 

V tx = 2(UyyV X - U X yVy) +UyV X y - U X Vyy (3-6) 
2 {Uyy U X ~\- 2 U X y H-y) 

is bi-Hamiltonian with respect to the standard Lie-Poisson structure on fl re g, together with \3. 5\) . 
where f = u y and a = v y . 

Proof. The simplest class of Casimir functions of this constant Poisson structure are linear 
combinations of the functionals f f dxdy and f udxdy. We will choose the Casimir function 



H o(f, a) = (a - f) dxdy. 

Js 1 xS 1 



The Hamiltonian vector field, Xh , with respect to the Lie-Poisson structure defines the following 
vector field 

ft = fx 

(3.7) 

at = 2f x + a x . 

Indeed, one obviously has (S a Ho, SfHo) = (—1, 1) (understood as an element of fl/R 2 ) and one 
then applies the definition (|2.4p . 

One thus looks for a function a) on g* eg such that its Hamiltonian vector field with 

respect to the constant Poisson structure satisfies 

X H X = X H > 

which leads to the following system of equation on the variational derivatives SfH\ and 5 U H\: 

-{SaH^ + ciSaH^ = f x 
2(8 a H 1 ) x -(5 f H 1 ) x + c(5 f H 1 ) y = 2f x + a x . 
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Introducing the first-order differential operator 

A = -d x + cd y , 

one shows by a simple straightforward calculation that following function: 

Hi(f, a) = [ (A-Hf x ) a + A"^/*) / - A~ 2 (f xx ) f) dxdy (3.8) 

is a solution of the above system. 

The Hamiltonian vector field Xh x is then as follows 

ft = A~Hfx)f x -A- 1 (fx X )f + c 2 A- 1 (f xy ) 

a t = A^if^^ + lA^if^a- A~ l {a xx )f -2A~ l (a x )f x 

-2 {A-\f xx ) - A- 2 (f xxx )) f - 4 {A-\f x ) - A~ 2 {f xx )) f x 
ci A- 1 {f xxxx )+c 2 (A- 1 (a xy ) + 2A- 1 (f xy )-2A~ 2 (f xxy )) 

In the same way as in [13], we substitute to this equation / = A(u) and a = A(v) and rewrite it 
in the following form: 

l^tx CUty — C {u X y 1l x ^y) ^2 l^xy 

VtX ~\~ CVty C (2 U XX Vy 2 U X y V X + U X V X y Uy V xx ) 

2 (u xx - A" (uxxx)) (u x -cuy) +4 (u x - A^iuxx)) (u xx - cu xy ) 
Cl u xxxx ~\~ c 2 (v X y ~\~ 2 U X y 2 A (u XX y)^) ■ 

It is very easy to check that, in the limit case c — ► oo, this system coincides with (|3.6p with 
exchanged notation for the variables (x,y) <-> (y,x). □ 

Theorem 13.41 implies the existence of an infinite series of first integrals in involution for the 
equation (jl.ip . as well as of an infinite hierarchy of commuting flows, see |16] . Section 5.5. 

Remark 3.5. 1) The special case c = in (|3.4p was considered in the details in [16]. This case 
is related to the equation (jl.2p . 

2) One can also choose a non-zero value of the first central charge c\ in (|3,4p . This will, 
however, only change the second equation in (|3.9p . 

3) Consider the first equation in (|3.9p . The term c 2 u xy can be removed by the transformation 
u i— ► u — x. Furthermore, the coordinate transformation (x, y) — ► (x,y + cx) leads to the 
following family: 

Ut x — C {Vj X y Uy l^yy ^x) ~\~ ^xx ^y ^xy ^x 

depending on c as parameter. This family gives one an interpolation between the equations (jl.ip 
and (jl.2p . but with zero central charge. 
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